Mori. Not. R. Astron. Soc. 000, |l||lO| (2012) Printed 22 April 2013 (MN M£X style file v2.2) 



Signatures of pulsars in the light curves of newly formed 
supernova remnants 
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ABSTRACT 

We explore the effect of pulsars, in particular those born with millisecond periods, 
on their surrounding supernova ejectas. While they spin down, fast-spinning pulsars 
release their tremendous rotational energy in the form of a relativistic magnetized 
wind that can affect the dynamics and luminosity of the supernova. We estimate the 
thermal and non thermal radiations expected from these specific objects, concentrating 
at times a few years after the onset of the explosion. We find that the bolometric light 
curves present a high luminosity plateau (that can reach lfj 43-44 erg/s) over a few 
years. An equally bright TeV gamma-ray emission, and a milder X-ray peak (of order 
10 40-42 erg/s) could also appear a few months to a few years after the explosion, as the 
pulsar wind nebula emerges, depending on the injection parameters. The observations 
of these signatures by following the emission of a large number of supernovae could 
have important implications for the understanding of core-collapse supernovae and 
reveal the nature of the remnant compact object. 

Key words: supernovae, superluminous supernovae, pulsars, pulsar winds, ultrahigh 
energy cosmic rays 



1 INTRODUCTION 

Core-collapse supernovae are triggered by the collapse and 
explosion of massive stars, and lead to the formation of black 
holes or neutron stars (see, e.g., Woosley et al. 20021. In 
particular, pulsars (highly magnetized, fast rotating neu- 
tron stars) are believed to be commonly produced in such 
events. The observed light curves of core-collapse supernovae 
present a wide variety of shapes, durations, and luminosi- 
ties, that many studies have endeavored to model, consider- 
ing the progenitor mass, explosion energy, radioactive nucle- 
osyntheis, and radiation transfer mechanisms in the ejecta 
(e.g.,|Hamuy|2003||Utrobin fc Chugai|2008| |Baklanov et al 
|2005||Kasen fc Woosley|2009| >. 

While they spin down, pulsars release their rotational 
energy in the form of a relativistic magnetized wind. The 
effects of a central pulsar on the early supernova dynamics 
and luminosity is usually neglected, as the energy supplied 
by the star is negligible compared to the explosion energy, 
for the bulk of their population. Some pioneering works have 
however sketched these effects Gaffet (1977a|bl; Pacini & 



Salvati ( 1973 1; Bandiera et al. ( 1984 1; Reynolds & Chevalier 



( 1984[), notably in the case of SN 1987A (|McCray et al.|1987 
Xu et al.|[T988l. Mo re recently, |Kasen fc Bildsten| ( j20To| 



Dessart et al. (12012) discussed that magnetars, a sub-class 



of pulsars born with extremely high dipole magnetic fields 
of order B ~ 10 14-15 G and millisecond spin periods, could 
deposit their rotational energy into the surrounding super- 
nova ejecta in a few days. This mechanism would brighten 
considerably the supernova, and could provide and expla- 
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nation to the observed superluminous supernovae (Quimby 
[2012] ). 

In this paper, we explore the effects of mildly magne- 
tized pulsars born with millisecond periods (such as the Crab 
pulsar at birth) on the light curves of the early supernovae 
ejecta. Such objects are expected to inject their tremendous 
rotational energy in the supernovae ejecta, but over longer 
times compared to magnetars (of order of a few years). In- 
deed, the spin-down and thus the timescale for rotational 
energy deposition is governed by the magnetization of the 
star. 

We estimate the thermal and non thermal radiations 
expected from these specific objects, concentrating at times 
of a few years after the onset of the explosion. We find that 
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the bolometric light curves present a high luminosity plateau 
(that can reach io 43-44 erg/s) over a few years, and that an 
equally bright TeV gamma-ray emission could also appear 
after a few months to a few years, from the acceleration 
of particles in the pulsar wind, depending on the injection 
parameters. A milder associated X-ray peak (of luminosity 
10 40-42 erg/s) could also be produced around the same time. 
The observations of these signatures by the following up of 
a large number of supernovae could have important implica- 
tions for the understanding of core-collapse supernovae and 
reveal the nature of the remnant compact object. 

These objects also present the ideal combination of pa- 
rameters for successful production of ultrahigh energy cos- 
mic rays (UHECRs, see |Blasi et al.|2000||Fang et al.|20"T2| ). 
The observation of such supernovae could thus be a fur- 
ther argument in favor of millisecond pulsars as sources of 
UHECRs, and a potential signature of an ongoing UHECR 
production. 

We first give, in Section [21 the list of quantities nec- 
essary for this analysis in the regimes of interest for the 
ejecta: optically thin or thick, and present a scheme of the 
early interaction between the pulsar wind and the supernova 
ejecta. In Section |3j we calculate the bolometric, thermal, 
and non thermal light curves of our peculiar supernovae. In 
Section [51 we briefly discuss available observations, and the 
implications for UHECR production. 
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Figure 1. Contour plot of the bolometric luminosity of supcr- 
nova+pulsar wind nebula systems at 1 yr after explosion (the 
fraction of wind energy converted into radiation, as defined in 
Section |3] is set to rj-y = 1), as a function of initial period P and 
magnetic field B. The various regimes for radiative emissions de- 
scribed in Section |2.1| are represented. The solid lines indicate 
pulsar spin-down timescale in seconds (Eq. [5j. The red dashed 
lines represent the pulsar population for which t p = tthim an( l 
separate naked and clothed pulsars (see text). The dotted lines 
represent t p = t^. 



2 SUPERNOVA EJECTA HOSTING A 

MILLISECOND PULSAR: PROPERTIES 

We note M e j and P c j as the mass and initial energy of the 
supernova ejecta. The pulsar has an inertial momentum 7, 
radius P» , initial rotation velocity Hi (corresponding initial 
period Pi = 27r/fii), and dipole magnetic field B. Numer- 
ical quantities are noted Q x = Q/10 x in cgs units, unless 
specified otherwise. 



2.1 Timescales 

For an ordinary core-collapse supernova, the ejecta expands 
into the circumstellar medium at a characteristic final ve- 
locity 



v e ] — «sn = 2 






1/2 



4.5 x10 s cm s^^M^ 



1/2 



(1) 
where M j j5 = M c j/5Mq. After a few expansion timescales 

t C x ~ R e j,i/vsN, where R e j,i is the radius of the star that 
led to the explosion, the ejecta enters into a stage of 
homologous expansion where its size scales as R — v c jt and 
its internal energy as E mt {t) ~ (E e j/2)(t ex /t). 

The ejecta is first optically thick to electron scattering. 
Noting k and p the opacity and density of the supernova 
envelope, one can estimate the optical depth of the ejecta: 
r = Rnp. Assuming a constant central supernova density 
profile (see |Matzner fc McKee||1999| and |Chevalier||2005| for 



more detailed modeling of the interior structure of super- 
novae) p = 3M c j/(4-kR 3 ), one can define the effective diffu- 



sion time (for thermal photons to cross the ejecta): 

1/2 



td = 



/ M c jK 

\4ttv c jC 



1.6 x 10 6 sMV 2 Kf 1/2 



ej.s-o.2 1 2xl0 9 cms-i 



-1/2 



(2) 
(3) 



with the opacity to electron scattering defined as 
K0.2 = /s;/(0.2g _1 cm ) for thermal photons. This sets 
the timescale of the supernova light curve, under the 
assumption that the opacity remains constant throughout 
the ejecta (no ionization effect), and in the absence of 
pulsar or 56 Ni heating. For more detailed computation of 



the these timescales, see, e.g., Kasen & Woosley (2009 



As the ejecta expands, it reaches a time ttMn when it 
becomes optically thin to electron scattering, for thermal 
photons (r = 1): 

1/2 

1.9 x 10 7 s ( - — 7 W^ ! t ) ■ (4) 



ithin 



3M G 



47rv 2 j 



2 x 10 9 cms" 



For the numerical estimates of v c j, we are using the final 
velocity of the ejecta after its modification by the shock at 
the interface between the pulsar wind and the initial ej ecta , 



for E, 



cj,51, M j i5 , 



and P; = 10" 3 s (see Eq.||]in Section [2~2 \ 



The pulsar spins down by electromagnetic energy losses 
that is transferred to the surrounding environment. The de- 
position of this energy happens over the spin-down timescale 
of the pulsar (Shapiro fc Teukolsky|l983): 



97c 3 



2B 2 RtQ,f 



3.1 x 10 7 s 74 5 Pr 3 2 Pr,6-Pi 2 - 



(5) 
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Figure 2. Scheme of the structures created by the interaction 
between the pulsar wind and the SN ejecta in the blast rest mass 
frame. 



We will consider two regimes for the calculation of ra- 
diative emissions from the ejecta: optically thin (t > tthin), 
and optically thick (i < tthin) for thermal photons. The de- 
position of pulsar rotational energy will have different effects 
on the supernova radiative emissions according to the opti- 
cal depth of the ejecta at time t p . Figure [l] pictures these 
various regimes. The red dashed lines represent the pulsar 



population for which t p 



on its left-hand side, most 



of the rotational energy of the pulsar is injected when the 
supernova ejecta is optically thin to electron scattering (the 
pulsar is naked). On the right-hand side of the red dashed 
line, the pulsar energy can enhance the luminosity of the 
supernova, as it is injected while the ejecta is still optically 
thick (the pulsar is clothed). 

2.2 Characteristics of the supernova ejecta and of 
the embedded pulsar wind nebula 

The interaction between the pulsar wind and the supernova 
ejecta leads to the formation of the following structures, il- 
lustrated in Fig. n\ a forward shock at the interface of the 
shocked and unshocked ejecta, and a reverse shock at the 
interface between the shocked and unshocked wind (com- 
monly called "termination shock"). The shocked material 
between the forward and the reverse shock constitutes the 
pulsar wind nebu la (PWN, e.g., |Ch evalier 1977 Chevalier 
|fc Fransson|1992| [Gaensler fc Slane||2006| ). 
The pulsar wind carries a total energy: 

and injects a luminosity (Shapiro & Teukolsky 1983) 



1.9 x 10 52 erg/ 4 5Pi 2 -3 , 



(6) 



Lp(t) 



E„ 



i p (1 + i/ip) 2 



(7) 



into the cold supernova ejecta. The evolution of the pul- 
sar luminosity over time, for magnetic dipole spin-down, is 
represented in Fig. 13] 

The characteristic velocity of the ejecta is not affected 
by the pulsar wind nebula expansion if E p <C E e j. How- 
ever, if the pulsar input energy overwhelms the initial ejecta 
energy E p S> E e j , the ejecta is swept up into the she ll at a 
final shell velocity Vi = (2E P /M ei ) 1/2 |Chevalier|2005| . Tak- 
ing into account these two extreme cases, one can estimate 
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Figure 3. Evolution of the pulsar luminosity L p as a function 
of time, for magnetic dipole spin-down. The pulsar has a dipole 
magnetic field of B = 10 13 G, and period Pi = 1, 10, 100 ms (in- 
creasing thickness). The vertical dashed lines correspond to the 
spin-down timescale t p for each intial spin period. 



the characteristic ejecta velocity as 

v C j =«sn(1 + £ p AE S n) 1/2 • ( 8 ) 

For E p <C -B e j , the evolution of the pulsar wind nebula 
takes place in the central part of the SN ejecta, where the 
density profile is nearly flat, with p oc t~ 3 (r/t)~ m . We will 
assume here that m = 0. For times t ^ t p where L p ~ 
Ep/tp, the radius of the pulsar wind nebula can then be 
expressed (Chevalier|1977l 



-Rpwn 



125 v^E p 

~99~ Mejtp 



1/5 



.6/5 



for t < t p , E p < E ei (9) 



Beyond the characteristic velocity usn, the density pro- 
file of the ejecta steepens considerably, reaching spectral in- 
dices b > 5 (e.g., |Matzner fc McKee||1999| |. For E p > E ej , 
the pulsar wind nebula expands past this inflection point 
and its size depends on whether the swept-up shell breaks up 
by Rayleigh- Taylor instabilities. Chevalier ( 2005 1 discusses 



that if the shell does not break up, the expansion is de- 
termined by the acceleration of a shell of fixed mass, thus, 
for t ^ t p , -E" p 3> E e j, and no shell disruption 



-RpWN 



_8__Ep_ 

15 M ej tp 



1/2 



,3/2 



2.2xl0 16 cm£j> e . 5 



1/2 



(10) 

fort = V (11) 



Otherwise, the evolution of the nebula is set by pressure 
equilibrium, and Rpwn oc t (6 ' b)/(5 ~ b) (for t < t p , E p > 
E G j). In the following, because the fate of the shell is unclear 



at this stage, we will use Eq. qlOb as an illustration. 

For t > tp, L p drops, and the swept-up material tends 
towards free expansion. One can roughly assume the relation 

-Rpwn(£ > tp) — -RpwN(tp)— , (12) 

t P 

where 7?pwN(t p ) is the size of the pulsar wind nebula in 
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time since explosion [days] 

Figure 4. Evolution in time of the radius, -RpwNi and the mag- 
netic field strength, Bp^Ni of a pulsar wind nebula, assuming no 
shell disruption (Eqs, [To] and |l2p and rye = 0.01, calculated for a 
SN ejecta with M B j = 5 Mq and E c j = 10 51 ergs -1 , embedding 



a pulsar with dipole magnetic field of B 
Pi = 1 ms. 



10 J G and period 



Eqs. ( 9|10 1 . More detailed modelings of the dynamical evo- 
lution of pulsar-driven supernova remnants can be found in 



Reynolds & Chevalier ( 1984 1 



The magnetic field strength in the pulsar wind nebula 
can then be estimated assuming a fraction of magnetization 
?7b of the luminosity injected by the wind (see Fig. H| 

Bpwn = (8717713 j L p {t')dt'j RpwN(ty' i/2 . (13) 

The value of r\B could vary between 0.01 — 1, according to 
pulsar wind nebulae. 



3 BOLOMETRIC RADIATION 

In what follows, we calculate the total radiation expected 
from the supernova ejecta+pulsar wind nebula. The evolu- 
tion of the ejecta is computed assuming a one zone core- 
collapse model. This approximation is debatable for times 
t < id, as the radiation should be mainly emitted in the 
central regions, close to the pulsar wind nebula, and not 
uniformly distributed as the matter over a single shell. This 
is not expected to be limiting for our study however, as we 
are most interested in the late-time light curves (a few years 
after the explosion), when the ejecta starts to become opti- 
cally thin. 

How much energy of the pulsar wind will be trans- 
formed into radiation depends on many factors such as the 
nature of the wind (leptonic, hadronic or Poynting flux 
dominated), the efficiency of particle acceleration and of 
radiative processes. In a first step, these conditions can be 
parametrized by setting a fraction rj-, of the wind energy 
E p that is converted to radiative energy (thermal or non 
thermal) in the pulsar wind nebula. 



Under the one zone model approximation, the radiation 
pressure dominates throughout the remnant, P = £i nt /3V, 
with V the volume of the ejecta. The internal energy then 
follows the law: 



il^ 



T)-,L v (t) — Lrad(i) 



(14) 



The radiated luminosity L ra d depends on the ejecta optical 
depth: 



£rad (t) 
4tt7? 2 

which yields 

irad(i) = 



-EintC 

(47t/3)i? 3 

E int C 

(47r/3)ri?3 



t>U 



t SC tthin 



3 -Bin 

Eintt 

' d 



t>u 



t Sj £thin 



(15) 
(16) 

(17) 
(18) 



where we note /3 e j = V e j/c. For t < tthin, we assumed that 
the totality of the luminosity 7? 7 L P deposited in the ejecta 
as photons is thermalized, and used the diffusion transport 



approximation Arnett (19801. In the optically thin regime, 



photons do not diffuse and propagate straight out of the 
ejecta. 



Equation ( 14 1 yields 

E in t(t) 



V~fE p 



1 + 3//3 oj 

- e d 
t 



hi 



t > tthi„(19) 



•/T P P 

(tp + x 



2 da; + _E cj i cx e ' 



t ^ tthin • 

The hypergeometric functions are noted: 

hi{x) = 2^(1,1 + 3/^,2 + 3/^,-3;) 
h 2 (x) = 2 Fi(2,l + 30 eh 2 + 3/0 eh -x). 



(20) 

(21) 
(22) 



Note that L ra d(t) ~ 7? 7 L p (t) for t > tthin- 

To calculate the total bolometric radiated luminos- 
ity, we add to L ra d(i) the contribution of the ordinary 
core-collapse supernova radiation Lsn(£)- Lstt(t) is cal- 
culated following Eq. (5) of Chatzopoulos et al. (2012), 
assuming an initial luminosity output of 10 42 erg/s, as is 
estimated by Woosley & Heger (2002) in their Eq. (41), for 



M c 



5 Mq and E c j — 10° ergs . Lsn only contributes 



when E p < Esn- 

Fig. [5] presents the bolometric luminosity radiated from 
the ejecta+PWN system for various sets of pulsar parame- 



ters. Again, the Arnett 1 1980[ ) approximation is not neces- 
sarily valid for t < td, where the radiation should not be 
distributed over the whole ejecta. Even with a rj-, < 10%, 
the plateau in the light curve a few years after the explosion 
is highly luminous, especially for P — 1 ms. This high lumi- 
nosity plateau stems from the injection of the bulk of the 
pulsar rotational energy a few years after the supernova ex- 
plosion. The luminosity is quickly suppressed for high B (for 
magnetar-type objects), due to the fast spin-down. Super- 
novae embedding isolated millisecond pulsars with standard 
magnetic field strengths would thus present unique radiative 
features observable a few years after their birth. 
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Figure 5. Evolution of the bolometric radiated luminosity of the supernova as a function of time. The pulsar has a dipole magnetic field 
of increasing strength as indicated, and period P; = 1 ms (left panel), and B = 10 13 G and increasing periods as indicated in the legend 
(right panel). The supernova ejecta has M j = 5 Mq and E c j = 10 erg. The gray lines give the evolution of the pulsar luminosity L p 
for each initial spin period. The gray dashed lines are the contribution of the ordinary core-collapse supernova to the radiated luminosity 
Lsn- We have assumed r; 7 = 1. The slight discontinuity is due to the numerical calculation of the integral in Eq. l |20| , and marks the 
transition between t < £ t hin an d t < tthin- 



Note however that the luminosity represented here is 
the bolometric one. The emission should shift from quasi- 
thermal to high energy after a few years, depending on the 
evolution of the opacity of the ejecta. The emission at dif- 
ferent energies is discussed in the next section. 



4 THERMAL/NON-THERMAL EMISSIONS 

The bolometric radiation calculated in the previous section 
stems from the re-processing of high energy radiation 
created at the base of the SN ejecta, in the PWN region. 
In the standard picture of PWN, high energy particles 
(leptons and hadrons) are injected at the interface between 
the pulsar wind and the ejecta, and radiate high energy 
photons (X rays and gamma rays). These high energy 
photons can be either thermalized if the medium (the 
PWN and/or the SN ejecta) is optically thick to these 
wavelengths, or can escape from the ejecta and be observed 
as a high energy emission, if the medium they have to 
propagate through is optically thin. In this section, we 
calculate in more detail the emission a few years after the ex- 
plosion, concentrating mainly on the case of a leptonic wind. 

Upstream of the termination shock, the energy of the 
pulsar wind is distributed between electrons and positrons, 
ions and magnetic fields. The fraction of energy imparted to 
particles is not certain, especially at these early times. Near 
the neutron star, the Poynting flux is likely to be the domi- 
nant component of the outflow energy. After many hundreds 
of years, observational evidence show that the energy repar- 
tition at the termination shock of pulsar wind nebulae is 
dominated by particles (e.g., Arons 2008). The conventional 
picture is thus that all but ~ 0.3 — 1% of the Poynting flux 



has already been converted into the plasma kinetic energy 
by the time the flow arrives the termination shock (Kennel 



fc Coroniti|1984a|b[[Emmering fc Chevalier|1 987; Bcg elman 



fc Li|1992 1, ^ 1% appearing to be a level required to repro 



duce the observed shape of the Crab Nebula ( Komissarov & 
Lyubarsky||2004 Del Zanna et al.pOOl l. How this transfer 



happens is subject to debate (see, e.g., Kirk et al.|2009 1 

Particles and the Poynting flux are injected in the pul- 
sar wind nebula at the termination shock. We will note 
the energy repartition between electrons and positrons, ions 
and the magnetic field in the pulsar wind nebula: L p — 
(rj c + rji+riB)L p . The ratio between rji and rj e is the subject 
of another debate (see e.g., Kirk et al.|2009[). However, var- 
ious authors (e.g. [Gelfand et al.|2009||Fang fc Zhang|2010 



Bucciantini et aL]|2011 Tanaka fc Takahara||20lT | seem to 
fit satisfactorily the observed emissions for various late time 
pulsar wind nebulae without adding any hadronic injection. 
We will thus focus on the emission produced for winds dom- 
inated by a leptonic component at the termination shock. 

Note that if protons are energetically dominant in the 



wind, Amato et al. 



( 2003 1 calculated that a large flux of 
neutrinos, gamma-rays and secondary pairs from p-p pion 
production should be expected from Crab-like pulsar wind 
nebulae around a few years after the supernova explosion. 
They estimate that the synchrotron emission from secon- 
daries will be negligible, while TeV photon and neutrino 
emission could be detectable by current instruments if such 
young objects were present in our Galaxy. 

Only 1% of the relativistic ions and magnetic fields com- 
ponents of the wind can be converted into thermal energy in 
the ejecta ( |Chevalier|| 1977| ) . This fraction can be amplified 
in presence of, e.g., Rayleigh- Taylor mixing, or high energy 
cosmic ray diffusion into the ejecta. 
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4.1 Pair injection in the PWN 



According to the original idea by Kennel & Coroniti ( 1984a I, 



the pair injection spectrum into the pulsar wind nebula 
should present a Maxwellian distribution due to the trans- 
formation of the bulk kinetic energy of the wind into ther- 
mal energy, and a non-thermal power-law tail formed by 
pairs accelerated at the shock. Hybrid and PIC simulations 



have shown indeed such a behavior (e.g., Bennett & Ellison 



1995|lDieckmann fc Bret|2009|[Sphkovsky|2008| >. |Spitkovsky' 
(2008) finds that 1% of the particles are present in this tail, 



with 10% of the total injected energy. The bulk of the par- 
ticle energy would then be concentrated around the kinetic 
energy upstream of the termination shock 



Cb 



kT c — 7 w m e c 
10 6 



5 x 10 11 eV 7 



(23) 
(24) 



with 7 W the Lorentz factor of the wind. The non-thermal 
tail would start around eb and continue up to higher ener- 
gies with a spectral index > 2. In practice, from a theoret- 
ical point of view, Lorentz factors as high as 7 W ~ 10 6 are 
difficult to reach, and current simulations are only capable 
of producing 7 W of order a few hundreds ( Spitkovsky||2008 

Sironi fc Spitkovsky|2009| . 

However, observationally, various authors (Kennel & 



Coroniti1|1984a| but also more recently, e.g., |Gelfand et al 
2009| |Fang fc Zhang|[2010l |Bucciantini et al.j|2011| |Tanaka 



fc Takahara|20lip demonstrated that the non-thermal radi- 
ation produced by the injection of either one single power- 
law or a broken power-law peaking around eb ~ 1 TeV, and 
extending up to PeV energies, could fit successfully the ob- 
served emission of various young pulsar wind nebulae. Such 
a high break energy implies either a high Lorentz factor for 
the wind 7 W ~ 10 5-6 , or an efficient acceleration mechanism 
enabling particles to reach 0.1 — 1 TeV energies. At higher 
energies, another acceleration mechanism has to be invoked 
to produce particles up to PeV energies. |Bucciantini et aT7| 
(2011 1 discuss that eb could possibly be viewed as a transi- 



tion energy between Type II and Type I Fermi acceleration 
from low at high energies. This would provide a physical 
explanation to the broken power-law shape, and alleviate 
the issue of the high wind Lorentz factor. At high energies, 
acceleration could also happen in the course of reconnec- 
tion of the striped magnetic field in the wind, at the ter- 
mination shock ( Lyubarsky|2003 Petri fc Lyubarsky|2 007). 
However, it is not clear yet whether this process can lead to a 
non-thermal particle distribution. One can expect additional 
particle acceleration in the wind itself, via surf-riding accel- 
eration (|Chen et al. 2002 Contopoulos & Kazanas 2002 



Arons 2002, 2003). This non-thermal component would not 



necessarily be processed when injected at the shocks if the 
particle Larmor radii are large compared to the size of the 
shock. 

In the following, we will assume that pairs are injected 
in the pulsar wind nebula following a broken power-law of 
the form 



diV 

~d7 



M) 



VeLp(t) 



(t/^)- a 


if £min ^ e < £b 


(e/e h y B 


if £b ^ £ ^ Emax 



(25) 

where a < 2 < /3, e m i n and e ma x are the minimum and maxi- 
mum cut-off energies respectively, and eb the peak of the in- 




thermal 
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phot + compt + pair of IC 
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1 10 100 1000 10000 

time since explosion [days] 

Figure 6. Timescales at play in the radiation emission of a PWN, 
for the same system as in Fig. [4] assuming an electron injec- 
tion break energy eb = 0.1 TeV, and t\b = 0.01. The dynami- 
cal timescale, id yn (black solid), and cooling timescales via syn- 
chrotron, 4 S yn (red solid), via self-compton, tic, syn (blue solid), 
and via IC off thermal photons (blue long dashed), are compared 
to thermalization timescales. The thermalization timescales via 
photoelectric absorption, Compton scattering, and pair produc- 
tion are in dashed lines (in red for the synchrotron photons and 
blue for the IC photons). Absorption by 7 — 7 interaction at high 
energies are in green lines (dot-dashed for the interaction between 
IC photons on synchrotron photons, and long dashed for the in- 
teraction of IC photons on the thermal photons of the SN ejecta). 
The gray shaded region corresponds to i < id, where the |Arnett| 
l|1980| approximation is not valid. And the dotted line indicates 
t = 'thin- 



jection distribution e(dA^/de) ~ 0.1 — 1 TeV. It is commonly 
assumed that eb oc 7 W oc y/L v {t), but such an assumption 
would imply very high wind Lorentz factors (> 10 9 ) at early 
times, that seem incompatible with the simulations and the- 
oretical models discussed above. It is likely that the Lorentz 
factor experiences a saturation above a certain value, and 
for simplicity, we will assume that eb is constant over time. 
For our purpose of deriving a rough estimate of the frac- 
tion of high energy emission that can escape the ejecta at 
early times, such an approximation will suffice. A thorough 
calculation of the emission spectrum would require time- 
dependent energy loss calculations for particles beyond the 
one zone approximation that we use here. |Del Zanna et a\\ 
( 2004 2006 1 have shown indeed that the high energy emis- 



sion is strongly affected by the details of the flow dynamics 
just downstream of the termination shock. 



4.2 Radiation by accelerated pairs 

The bulk of the electron distribution will predominantly ra- 
diate in synchrotron and experience inverse Compton (IC) 
scattering off the produced synchrotron photons. The cool- 
ing timescales of these processes, as well as the dynami- 
cal timescale tdyn = -Rpwn/c of the PWN are indicated 
in Fig. [6] Inverse Compton scattering off the thermal pho- 
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Figure 7. Evolution in time of radiative quantities of the same system as in Fig. El with T]b = 0.01, r) e = 1 — rjg, and for a particle 
injection break energy of £b = 0.1 TeV (left) and eb = 1 TeV (right). Blue dashed line: characteristic energy of synchrotron radiation; red 
dashed lines: characteristic energy of photons produced via self-comptonization of synchrotron photons. Thin black solid line: luminosity 
injected by the pulsar. Red and blue solid lines: luminosity emitted in synchrotron and IC radiations. Dotted blue line: optical depth 
of ejecta of composition 60% H, 30% He and 10% C, to photons of energy hu c . Dotted red lines: optical depth of ejecta to photons of 
energy v lc . 



tons of the ejecta and off the Cosmic Microwave Background 
(CMB) are negligible compared to the former two processes. 
Figure [6] also demonstrates that the cooling timescale of IC 
scattering off the thermal photons of the SN ejecta is much 
longer than the timescale for self-comptonization of the syn- 
chrotron emission. This estimate includes only the contri- 
bution of the thermal photons of the standard supernova 
ejecta, as the thermalization of the non-thermal components 
described here happen on larger timescales, in the optically 
thin regime which is of interest to us. 

The synchrotron cooling timescale of accelerated elec- 
tron reads 



Csvn — 



4CTTCct/fl 



(26) 



with Ub = -Bpwn/871"- At the early stages that we con- 
sider here, the characteristic energy of radiating particles 
is e c (t) = £b- The characteristic synchrotron radiation fre- 
quency can be expressed 



u c (t) = 0.29 



3eB F 

47rm3c 5 



>(t). 



(27) 



Accelerated electrons also scatter off these synchrotron pho- 
tons by IC, producing photons at energy 



I'lC 



"KN 



Jell 
if V c < l^KN 



if V c > !/kn , 



with a cooling timescale 

ilC.syn = 



3mc 4 



AlJTCtcUs- 



(28) 

(29) 

(30) 



where U syn is the synchrotron photon energy density. Elec- 
trons radiate in synchrotron and self-compton processes 



with the following power ratio: Pic/P S \ 



U sy n/U B - As- 



suming that the energy of the accelerated electron popula- 
tion is concentrated in its peak energy eb, this implies syn- 
chrotron and IC luminosities of L syn = r)Br] e /(r]B + »7e)L p 
and Lie = Ve/(VB+Ve)^p respectively. Obviously, the value 
of r\B has an impact on the synchrotron emission, but not 
on the IC emission. 

Figure W\ presents the evolution in time of the luminosi- 
ties L p , Lsyn, and Lie, as well as the emission frequencies 
v c and u\c- The IC radiation is mostly emitted at the break 
energy of the injection of electrons. The synchrotron emis- 
sion spans from gamma/X-ray (until a few years) to optical 
wavelengths (after thousands of years). At the time of inter- 
est in this study, X-rays are thus mainly emitted between 
0.1 - 100 keV for e b = 0.1 TeV, and around 100 keV - 1 GeV 
for e b = 1 TeV. 



4.3 Thermalization in the ejecta 

The X-ray opacity (~ 0.1 — 100 keV) is dominated by pho- 
toelectric absorption in metals. Above ~ 100 keV, very hard 
X-rays and gamma-rays experience predominantly Compton 
scattering, and pair production above ~ 10 MeV. The opac- 
ities of these processes for various atomic media are given in 
Fig. [8] At a given time t, the optical depth of the ejecta to 
the characteristic synchrotron photon emission v c reads 



yn(t) = V e jtKv c (t)p(t) 



(31) 



The dominant thermalization process for the TeV IC radi- 
ation is pair production by 7 — 7 interactions (see Fig. p|. 
The timescale for thermalization via this process is however 
slightly longer than the dynamical time; hence the ejecta 
appears mostly optically thin to this radiation (see Fig. uv. 
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Figure 8. Atomic scattering opacities of high energy photons 
on H, He, C, and a mixture of composition mimicking that 
of a type II core-collapse supernova ejecta (60% H, 30% He 
and 10% C). The black dashed line indicates the contribution 
of Compton scattering in the latter composition case. From 
http://henke.lbl.gov/optical_constants/ for 30 eV ^ hv < lkcV 
and http://physics.nist.gov/cgi-bin/Xcom/ for 1 keV ^ hv < 
100 GeV. 



We note tic f° r the optical depth of the ejecta to the IC 
emission. 

The luminosity in the characteristic energies hu c and 
hv\c after thermalization, noted respectively Lx(t) and 
L~ t (t), and the luminosity in thermal photons, L t h(t), are 
calculated as follows 

L x (t) = L Byn (v c ,t)e- T ^ (t) (32) 

L 7 (£) = L ic (v ic ,t)e- Tic{t) (33) 

Lth(t) = L TaA (t) - L x (t) - L~,(t) . (34) 

Figure [9] presents the thermal emission (black), X-ray 
emission (blue dotted) at hv c ~ 0.1 — 100 keV for eb = 
0.1 TeV (left panel) and ~ 100 keV - 1 GeV for e b = 1 TeV 
(right panel), and 0.1 — 1 TeV gamma ray emission (red 
dashed) expected from a SN ejecta with M e j = 5 Mq and 
E e j — 10 51 ergs -1 , embedding a pulsar with dipole magnetic 
field of B — 10 13 G and period Pi = 1,3,10 ms (increasing 
thickness), assuming 77s = 0.01, r/ e = 1 — r\s, and a break 
energy e b = 0.1 TeV (left) and e b = 1 TeV (right). 

A decrease in flux is expected in the thermal com- 
ponent after a few months to years, when the ejecta be- 
comes optically thin to gamma-rays. For a low break energy 
(eb = 0.1 TeV) the thermal component can then recover, as 
the X-ray emission vanishes, because of the increase of the 
ejecta optical depth for lower energy photons. One robust 
result is that, in both break energy cases, for fast pulsar 
rotation periods Pi ^ 3 ms, the associated gamma-ray flux 
around 0.1 — 1 TeV emerges at a level that should be de- 
tectable at a few tens of Mpc, and remains strong over many 
years. 



5 DISCUSSION, CONCLUSION 

We have estimated the thermal and non thermal radiations 
expected from supernova ejectas embedding pulsars born 
with millisecond periods, concentrating at times a few years 
after the onset of the explosion. The bolometric light curves 
should present a high luminosity plateau (that can reach 
> 10 43 erg/s) over a few years. A more detailed emission cal- 
culation considering the acceleration of leptons in the pulsar 
wind nebula region shows that an X-ray and a particularly 
bright TeV gamma-ray emission (of magnitude comparable 
to the thermal peak) should appear around one year af- 
ter the explosion. This non thermal emission would indicate 
the emergence of the pulsar wind nebula from the supernova 
ejecta. 

The light curves calculated in this paper are simple es- 
timates, that do not take into account second order effects 
of radioactive decay of 56 Ni, recombination, etc. (see, e.g., 
Kasen & Wooslcy 2009). The non-thermal components are 
also evaluated assuming that all the leptonic energy is con- 
centrated in one energy bin. A more detailed analysis should 
be conducted, taking into account the shape of the spectra 
and its evolution in time, in order to get a more accurate rep- 
resentation of the emission, and for a thorough comparison 
with observational data. Depending on the spectral indices, 
a non monoenergetic electron injection spectrum could lead 
to a decrease of the peak luminosity of one order of magni- 
tude. 

Our computation of the evolution of the PWN (radius, 
magnetic field) is also basic, and could benefit from more 
thorough estimations. Our toy model suffices however in the 
scope of this study, where the aim is to demonstrate the im- 
portance of multi-wavelength follow ups of SN lightcurves. 
We also assumed a relatively high magnetization r\B of the 
wind at the termination shock, following estimates that re- 



produce the features of the Crab nebula (Komissarov & 
Lyubarsky|2004| |Del Zanna et al.|2004 l. 



Several earlier works treat some of the aspects invoked 
above in more detail. For example, in the context of evolu- 
tion of pulsar wind nebulae, early works by [Pacini & Salvati 



|T973l);[Ree"s fc Gunn|(|l974|);|Bandiera et al 



& Panagia ( 1980 ) ; Reynolds & Chevalier ( 1984 1 take into ac- 



1984) 



Weiler 



count the detailed evolution of particle energy distribution 
and radiation spectrum. Most of these works aim at calcu- 
lating radiation features of observed plerions, a few hundred 
of years after the explosion. However, their modeling at ear- 
lier times, especially in the work of Reynolds & Chevalier 
(19841, lays the ground for the more detailed calculations 



that should be performed in our framework. 

The level of synchrotron emission predicted here can 
thus be viewed as optimistic values. However, the gamma- 
ray flux that is predicted does not depend on the magnetiza- 
tion, and remains fairly robust to most parameter changes. 

Currently, only a handful of supernovae have been fol- 
lowed over a period longer than a year, and no object, except 
for SN 1987A, has been examined in X-rays or TeV gamma- 
rays a year after the explosion. Among the objects that have 
been followed in optical bands, SN 2003ma ( Rest et al.| 2011) 
has an abnormally bright luminosity at the peak, and a long 
bright tail over many years. The six type II supernovae fol- 
lowed by Otsuka et al. ( 2012 I present various shapes of light 
curves, and a cut-off in the thermal emission after a few 



© 2012 RAS, MNRAS 000,[T]fIO 



Signatures of pulsars in the light curves of newly formed supernova remnants 9 



0.1 TeV 








e b =1TeV 




1(J |J 










10 44 






r ****^ X-rays x10 

A "** *— ^ gommo-rays 

T~\ " ""~ - _ 




3 


1Q 43 


f ^\ '"'•■ 


o 
J 


10 42 


K \. ~~ 




10 4 ' 


v ^-^^____^_ 




i' /X 




m 40 


1; 1 \. 

i : , i , i TS. , i , , , i , , , i , , . i , , , i . , , i , 



200 400 600 800 1000 1200 1400 
t [days] 



200 400 600 800 1000 1200 1400 
t [days] 



Figure 9. Thermal emission (black solid lines), non thermal X-ray emission (blue dotted) at energy indicated in Fig. W\ (emitted mostly 
between 0.1 — 100 keV for ej, = O.lTeV, and around 100 keV— 1 GeV for t\, = ITeV), and non thermal gamma-ray emission (red dashed) 
at 0.1 TeV (left) and 1 TeV (right), expected from a SN ejecta with M e j = 5 Mq and E e j = 10 51 erg s — 1 , embedding a pulsar with dipole 
magnetic field of B = 10 13 G and period P; = 1,3, 10ms (increasing thickness), assuming r/g = 0.01, rj e = 1 — r}g, and a break energy 
eb = 0.1 TeV (left) and ej, = ITeV (right). Caution: for visibility, the X-ray luminosity is multiplied by 10 on the left panel, and by 10 
on the right panel. The slight discontinuities are numerical artifacts. 



years. Our study demonstrates that the features in these 
light curves could also be due to the energy injection from 
the pulsar, that could compete with the other processes that 
are more commonly considered, such as the light echo of the 
peak luminosity, or the radioactive decay of 56 Ni. An associ- 
ated X-ray and TeV gamma-ray emission emerging around 
a few months to a year after the explosion would constitute 
a clear signature of pulsar rotational energy injection. It is 
also interesting to note that the emergence of a pulsar wind 
nebula has been recently reported from radio observations 



for SN 1986J (Bietenholz et al. 20101, though over longer 



timescales than predicted for the objects studied in this pa- 
per. 



Some authors (Katz et al. 2011 Svirski et al. 20121 



have discussed that shock breakouts from stars surrounded 
by a thick wind could lead to bright X-ray peak after a few 
months, similar to the signal discussed in this paper. This de- 
generacy can be overcome by the observation of the gamma- 
ray signal, which should be absent in the shock breakout 
scenario. Detailed analysis of the respective X-ray spectra 
should also help distinguish the two scenarios. 

The follow up of bright type II supernovae over a few 
years after the explosion in different wave bands would thus 
reveal crucial information on the nature of the compact rem- 
nant. These suveys should be made possible with the advent 
of optical instruments such as LSST, and the use of powerful 
instruments for transient event detection, such as the Palo- 
mar Transient Factory or Pan- STARR. The bright X-ray 
signal should be detected by NuSTAR for supternovae out 
to redshifts z ~ 0.5, and the even brighter gamma-ray signal 
could be observed by HESS2, by the future Cerenkov Tele- 
scope Array (CTA), and by HAWK which will be the choice 
instrument to explore the transient sky at these energies. 
For CTA, an adequate survey of the sky outside the Galactic 



plane could spot gamma-ray sources of luminosity 10 43 erg/s 
as predicted by this work within a radius of ~ 150 Mpc (G. 
Dubus, private comm.). Assuming a gaussian pulsar period 
distribution centered around 300 ms as in |Faucher-Giguere| 
& Kaspi (2006b, implies that 0.3% of the total population 



has spin periods < 6 ms. With this estimate, one could find 
4 bright sources within 150 Mpc. This is consistent with the 



numbers quoted in early works by Srinivasan et al. ( 1984 1 



|Bhattacharya| ( |1990[ ) . These authors estimated the birthrate 
of Crab- like pulsar-driven supernova remnants to be of order 
1 per 240 years in our Galaxy. 

Pulsars born with millisecond periods embedded in 
standard core-collapse supernova ejectas, as described in this 
paper, are promising candidate sources for ultrahigh energy 
cosmic rays ( |Fang et al.||2012} |2013[ ). In the framework of 
UHECRs, an injection of order 1% of the Goldreich- Julian 
density into ions would suffice to account for the observed 
flux, assuming that 1% of Type II supernovae give birth to 
pulsars with the right characteristics to produce UHECRs 
(i.e., pulsars born wi th millisecond pe riods and magnetic 
fields B ~ 10 12 ~ 13 G, rFang et al.|2012| ). The observation of 
the peculiar light curves predicted here could thus provide 
a signature for the production of UHECRs in these objects. 
Though no spatial correlation between arrival directions of 
UHECRs and these supernovae is expected (because of time 
delays induced by deflections in magnetic fields), an indica- 
tion of the birth rate of these supernovae could already give 
direct constraints on this source model. Photo-disintegration 
and spallation of accelerated nuclei in the supernova ejecta 
could also lead to an abundant high energy neutrino produc- 
tion ( Murase et al.|200"9 consider such a neutrino production 
in the case of magnetars instead of fast-rotating pulsars), 
that could be detected with IceCube, and correlated with 
the position of identified peculiar supernovae. 
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